Quantum Monte Carlo calculations of excited states in A = 6 — 8 nuclei 
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A variational Monte Carlo method is used to generate sets of orthogonal trial functions, 3 , t(J 71 ,T), 
for given quantum numbers in various light p-shell nuclei. These are then used as input to Green's 
function Monte Carlo calculations of first, second, and higher excited (J T , T) states. Realistic two- 
and three-nucleon interactions are used. We find that if the physical excited state is reasonably 
narrow, the GFMC energy converges to a stable result. With the combined Argonne vis two- 
nucleon and Illinois-2 three-nucleon interactions, the results for many second and higher states in 
A — 6-8 nuclei are close to the experimental values. 

PACS numbers: PACS numbers: 21.10.-k, 21.45.+V, 21.60.Ka 



I. INTRODUCTION 

Quantum Monte Carlo methods have proved to be very 
accurate and powerful tools for studying the structure of 
light nuclei with realistic two- and three-nucleon inter- 
actions. In a series of papers, we have calculated about 
40 ground and low-lying excited state energies of differ- 
ent (J* ,T) quantum numbers in A = 6-10 nuclei with 
an accuracy of ~ 1-2% E H 0, . The first step is 
a variational Monte Carlo (VMC) calculation to find an 
optimal trial function, ^T(J n ,T) for a given state. The 
'Ft, which is antisymmetric by explicit construction, is 
then used as input to a Green's function Monte Carlo 
(GFMC) algorithm, which projects out the lowest-energy 
eigenstate by a propagation in imaginary time, r. The al- 
gorithm preserves the quantum numbers of although 
it may introduce symmetric noise into the propagated 
'F(t). In principle, ^B(t) for large r approaches the ex- 
act lowest-energy eigenstate with the specified quantum 
numbers. 

VMC calculations of second and higher excited ( J*", T) 
states have also been made, because a major step in the 
preparation of an optimal "Fy is a diagonalization in the 
small single-particle basis of different possible spatial- 
symmetry states 0, S ■ GFMC calculations of second 
or higher (J 7r ,T) states have not been attempted previ- 
ously in nuclear physics, under the expectation that any 
small contamination of the excited trial state by the true 
first state would drive the calculated energy below its 
proper value. However, as is described in this paper, we 
have found that it is possible to obtain useful GFMC pre- 
dictions of multiple states with the same quantum num- 
bers. 

In this paper we report such GFMC calculations of 
second and higher (J* ,T) states in light (A = 6-8) p- 
shell nuclei. VMC calculations of these states were re- 
ported previously 0, HI, but we have made significant 
improvements in the trial functions since then. Details 



of the recent VMC work are discussed in Sec. II. For the 
GFMC calculations, we find that the propagated energy, 
E(t), for many states is stable and a useful excitation 
energy can be extracted. GFMC calculations for ground 
states are reviewed in Sec. Ill, where we also describe an 
improved algorithm for propagation with three-nucleon 
interactions that is both faster and more accurate than 
what we had used previously. The GFMC algorithm for 
higher excited states is described in Sec. IV, along with 
various orthogonality tests. Numerical results are given 
in Sec. V for the Argonne vig (AV18) and simplified v§' 
(AV8') two-nucleon (NN) interactions fi], and for AV18 
with either the Urbana IX (UIX) or Illinois-2 (IL2) three- 
nucleon (NNN) interactions Q, Q added. We find, con- 
sistent with our studies of first (J 71 , T) states, that the 
AV18/IL2 Hamiltonian gives a good representation of the 
experimental spectrum; the RMS deviation from 36 ex- 
perimental energies with 6 < A < 8 is only 0.60 MeV. 
Some concluding remarks are given in Sec. VI. 



II. VMC TRIAL FUNCTIONS 



The VMC trial function, * T (J Tr ;T), for a given nu- 
cleus, is constructed from products of two- and three- 
body correlation operators acting on an antisymmet- 
ric single-particle state with the appropriate quantum 
numbers. The correlation operators reflect the influ- 
ence of the interactions at short distances, while appro- 
priate boundary conditions are imposed at long range. 
The contains variational parameters that are ad- 
justed to minimize the energy expectation value, Et — 
(^t\H\^ t )/(^t\^t), which is evaluated by Metropolis 
Monte Carlo integration. 
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A good variational trial function has the form 



|*t> = 



i<j<k 



mi 



Ui. 



I*/) - (1) 



The Uij and U™ 1 are non-commuting two- and three- 
nucleon correlation operators, S indicates a symmetric 
sum over all possible orderings, and ^ j is the fully anti- 
symmetric Jastrow wave function. For s-shell nuclei the 
latter has the simple form 



i<j<k 



i<3 



\$ A (JMTT 3 )} . (2) 



Here fc^ij) an d f^ k are central two- and three-body 
correlation functions and <$>a is a Slater determinant in 
spin-isospin space, e.g., for the a-particle, $4(0000)) = 

The correlation operator Uij includes spin, isospin, and 
tensor terms: 



u a = X! u p( r v)°ij - 

p=2,6 



(3) 



where the 0\~ ' — [1, avcry, Sy 55 [1, U-Tj] are the same 
static operators that appear in the AW potential. The 
functions f c (r) and u p (r) are generated by solving a set 
of six coupled differential equations with embedded vari- 
ational parameters 0: two single-channel equations in 
1 S and 1 P waves, and two coupled-channel equations in 
3 S and 3 P waves. The U™ 1 has the spin-isospin struc- 
ture of the dominant parts of the NNN interaction as 
suggested by perturbation theory. 

For p-shell nuclei, the Jastrow wave function, Vl'j, is 
necessarily more complicated. It starts with a sum over 
independent-particle terms, <5> A , that have 4 nucleons in 
an a-like core and (A — 4) nucleons in p-shell orbitals. 
We use LS coupling to obtain the desired JM value of a 
given state, as suggested by standard shell-model stud- 
ies §|. We also need to specify the spatial symmetry [n] 
of the angular momentum coupling of the p-shell nucle- 
ons |9j. Different possible L5[n] combinations lead to 
multiple components in the Jastrow wave function. This 
independent-particle basis is again acted on by products 
of central pair and triplet correlation functions, but now 
they depend upon the shells (s or p) occupied by the 
particles and on the LS[n] coupling: 



i*./) = A [ n /&*][ n /"M e (a*[n][ n /. 

[ i<j<k i<j<4 LS[n] k<4<l<A 

x [ II fp P Sln] ( r im)] \® A (LS[n]JMTT 3 ) 12U:5 ... A )} \ 

4<l<m<A \ 



Km) 



(4) 



The operator A indicates an antisymmetric sum over all possible partitions of the A particles into 4 s-shell and (A — 4) 
p-shell ones. The pair correlation for both particles within the s-shell, f ss , is similar to the f c of the a-particle. The 
pair correlations for both particles in the p-shell, fp P S ^ n \ and for mixed pairs, fs P S ^ n \ are similar to / ss at short 
distance, but their long-range structure is adjusted to give appropriate clustering behavior, and they may vary with 
LS[n}. 

The single-particle wave functions $a are given by: 



\$A(LS[n]JMTT 3 ) 123 4: 5 ...A} 



|<I>,(()n()m l234 [ J] cb^- 

4<1<A 



<Pi s[n] (R«i) 



(5) 



4<1<A 



L Ml [n] 



n ^\ ms ) 



4<1<A 



SM S 



[ n 



JM 



4<1<A 



TT 3 



The 4> p (R a i) are p-wave solutions of a particle in 
an effective a — N potential that has Woods-Saxon and 
Coulomb parts. They are functions of the distance be- 
tween the center of mass of the a core and nucleon I, and 
may vary with LS[n]. The depth, width, and surface 
thickness of the single-particle potential are additional 
variational parameters of the trial function. The overall 



wave function is translationally invariant, so there is no 
spurious center of mass motion. 

The 0LS{n] mixing parameters of Eq. are deter- 
mined by a diagonalization procedure, in which matrix 
elements 

E T ,ij = <*t (ft) \H\* T (ft)) , (6) 
N T ,a = (*t(A)|*t(/?j)} , (7) 
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are evaluated using trial functions ^riPi) = ^t(A = 
1,/3,-^i = 0). Although the <P A (LS[n]JMTT 3 ) are or- 
thogonal due to spatial symmetry, the pair and triplet 
correlations in 4"^ make the different LS'[n] components 
nonorthogonal, so a generalized eigenvalue routine is nec- 
essary to carry out the diagonalization. 

TABLE I: f3 L s[ n ] components for 6 He states, listed in order 
of increasing excitation. (L = S or D as appropriate.) 

(jyr) x l[2] 3 p[ii] 

(0+;l) 0.974 -0.228 

(2+;l) 0.922 0.386 

(2+;l) -0.388 0.922 

(1+;1) 1. 

(0+;l) 0.232 0.973 

TABLE II: PlsM components in 6 Li states, listed in order of 
increasing excitation. 

(J";T) 3 S[2] 3 D[2] ^[11] 

(l + ,0) 0.986 0.138 0.098 

(3 + ;0) 1. 

(2 + ;0) 1. 

(l + ,0) -0.109 0.964 -0.242 

(l + ,0) -0.134 0.231 0.964 



TABLE V: Plsm components for He states, listed in order 
of increasing excitation. (L = S or D as appropriate.) 

{J W ;T) *L[22] 3 P[211] 
(0 + ;2) 0.794 -0.608 
(2 + ;2) 0.928 -0.373 
(l+;2) 1. 
(0 + ;2) 0.610 0.792 
(2+;2) 0.377 0.926 

are included. The added components have little effect 
on the energies of the lowest states of given ( J 7r ,T), but 
their inclusion can be important for higher excited states; 
they are also expected to be important for some elec- 
troweak transitions. The diagonalizations were done for 
the AV18/UIX Hamiltonian and the resulting Plsm used 
for all the Hamiltonians reported here; a few tests have 
shown insignificant dependence of the /3ls[u] on Vyfc. 

Additional improvements in the trial wave have been 
made by refining the detailed shapes of the u p (r), f sp {r), 
and f pp (r) functions; these improvements are described 
in more detail with the numerical results in Sec.V. 



III. GFMC METHODS FOR GROUND STATES 



TABLE III: Pls[u] components for 
of increasing excitation. 

(J^T) 2 P[21] 2 D[21] 



He states, listed in order 
4 S[111] 



(f 

V 2 
^ 2 
1 2 



2 ) 
' 2> 

2 ) 
■ 21 

' 2> 

2 ) 

' 2> 

2) 
• 2 > 



0.837 
1. 

-0.392 
0.373 



0.515 -0.182 



1. 

0.805 
-0.317 



0.445 
0.872 



TABLE IV: Pz,s\„] components for 7 Li states, listed in order 

^ or F as appropriate.) 





2 L[3] 


4 P[21] 


4 D[21] 


2 P[21] 


2 D[21] 


2 S[111] 


(3- 1) 

\2 ' 21 


0.995 


0.086 


0.026 


0.005 


-0.050 




fl" I) 

\2 ' 21 


0.988 


-0.003 


-0.098 


-0.116 




-0.024 


(1- 1) 

\2 ' 21 


0.990 




0.138 








(5- 1) 

\2 ' 2> 


0.988 


0.114 


0.072 




-0.079 




(I" 1) 

V2 ' 2> 


-0.079 


0.957 


-0.058 




0.274 




V2 ' 2> 


-0.084 


0.952 


0.250 


-0.087 


0.132 




V2 ' 2> 


0.049 


0.841 


-0.111 


0.523 




-0.069 


(1- 1) 

V2 ' 2> 


-0.129 




0.992 








V2 ' 2> 


-0.088 


-0.009 


0.968 




0.235 




(5- I) 
\2 ' 2> 


0.137 


-0.263 


-0.224 




0.928 





Green's function Monte Carlo calculations of light nu- 
clei have previously been performed for ground states 
of light nuclei 0, 0, 0, 0, E| and for the lowest-energy 
states with distinct quantum numbers. In this section 
we briefly review the GFMC method as applied to light 
nuclei; Refs. 010 should be consulted for detailed de- 
scriptions of the method and tests of its reliability. We 
then report new algorithmic developments which signifi- 
cantly increase the speed of these calculations. These im- 
provements result in a much more efficient treatment of 
NNN interactions, particularly the very small, but com- 
putationally expensive, terms when three nucleons are 
well separated. In the next section we discuss the exten- 
sions required to treat excited states of the same quantum 
numbers. 

Starting with the trial function GFMC provides a 
means of computing a wave function propagated in imag- 
inary time: 

*(r) = exp[-(H - E )t]^> t = exp[-(H - E )AT} n ^ T , 

(8) 

by stringing together a series of small time ( At = t/N) 
evolution operators. The energy, E(t), is evaluated using 
a "mixed" expectation value between $>t and \&(r): 



The results for the (3i,s \n\ c oefficients are shown in Ta- 
bles [IHVTTI Tables |l| |ffl| and |V| show the complete 
set of possible p-shell states for 6 He, 6 Li, 7 He, and 8 He, 
respectively. Tables IIVI I VII and IVIII show the lowest- 
lying p-shell states for 7 Li, 8 Li, and 8 Be. In previous 
work only the two most spatially symmetric LS[n] 
components were included in 8 Li and 8 Be, but now all 



E(t) = <#(T)) M ixed 



(^(t)\H\^ t ) 
(*(t)|¥t) 



(9) 



In principle, E(t) is an upper bound to the energy of the 
lowest eigenstate that is not orthogonal to ^>t- 

The ^(t) is represented by a vector function of R, 
where the components of the vector are the amplitudes 
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TABLE VI: /3ls[„] components for 8 Li states, listed in order of increasing excitation. (L = P, D or F as appropriate.) 



(J n 


T) 


3 P[31] 


3 D[31] 


3 F[31] 


: L[31] 


3 S[22] 


3 D[22] 5 P[211] 


3 P[211] 


^[211] 


(2+ 


1) 


0.948 


-0.248 


-0.094 


0.122 




-0.116 


0.047 


0.019 




(1 + 


1) 


0.772 


-0.299 




0.524 


-0.091 


-0.074 


0.120 


-0.112 


-0.005 


(3+ 


1) 




0.947 


0.126 


-0.257 




-0.061 


-0.135 






(0+ 


1) 


0.995 














-0.100 




(1 + 


1) 


-0.588 


-0.302 




0.733 


0.123 


0.089 


-0.025 


-0.031 


-0.009 


(2+ 


1) 


0.281 


0.908 


0.297 


-0.038 




0.040 


-0.056 


0.057 




(2+ 


1) 


-0.079 


-0.128 


0.533 


0.823 




-0.092 


-0.053 


0.063 




(1 + 


1) 


0.051 


0.873 




0.412 


-0.221 


0.065 


-0.045 


-0.064 


0.078 


(3+ 


1) 




-0.153 


0.945 


-0.125 




-0.229 


0.126 






(4+ 


1) 






1. 















TABLE VII: /3i,s[ji] components for 8 Be states, listed in order of increasing excitation. (L = S, D or G as appropriate.) 





T) 


l L[4] 


3 P[31] 


3 D[31] 


3 F[31] 


5 S[22] 5 D[22] 


*L[22] 


3 P[211] 


(0+ 


0) 


0.998 


-0.055 








-0.034 


0.016 


0.015 


(2 + 


0) 


0.998 


0.026 


0.041 


-0.020 


0.002 


0.029 


0.007 


0.008 


(4+ 


0) 


0.999 






0.052 




-0.015 






(2 + 


0) 


-0.011 


0.931 


-0.169 


0.087 


0.284 


-0.062 


-0.101 


-0.055 


(1 + 


0) 




0.957 


0.237 






-0.165 




0.018 


(1 + 


0) 




-0.249 


0.959 






0.098 




0.097 


(3 + 


0) 






0.936 


0.287 




0.207 






(4 + 


0) 


-0.018 






0.691 




0.723 






(2 + 


0) 


-0.038 


-0.006 


0.825 


0.230 


0.499 


0.050 


0.111 


0.039 


(0+ 


0) 


0.029 


0.926 








-0.311 


0.181 


-0.105 


(3+ 


0) 






-0.127 


0.782 




-0.610 






(2+ 


0) 


0.016 


-0.301 


-0.376 


-0.257 


0.782 


0.148 


-0.011 


-0.261 



of each of the distinct spin-isospin states in the system. 
The propagator 

G Q/3 (R,R'; At) = (R, a\ exp[— (i? - Eq)At] |R', /3) 

(10) 

depends upon the 6 A spatial coordinates R and R', as 
well as the spin-isospin states a and f3. It is calculated 
with leading errors of (At) 3 . The full wave function is 
then obtained as a product over many small time steps: 

*(R»,t) = y , G(R n ,R„_ 1 )---G(R 1 ,R ) 



* r (Ro)dRr, 



• dRo, 



(11) 



where we have omitted the spin-isospin labels for sim- 
plicity. 

A central ingredient in a GFMC calculation is the 
propagator G(R, R'; At). It is desirable to use as large a 
time step At as possible to speed up the calculation, how- 
ever the maximum time step is limited by the accuracy 
with which we can evaluate the propagator G(R, R'; At). 
The short-time propagator typically used for NN inter- 
actions is 

G a)/9 (R,R';AT) - 



G (R,R')(a| 



9ij { r ij i r ij ) 
90,ij { r ij > v ' ij ) 



(12) 



where the free-particle propagator Go and the free two- 
nucleon propagator go.ij are simple Gaussians and 
the interacting pair propagator, 



; At) = (ry|exp[-AT(r<j 



,), (13) 



is easily calculable from the two-body potential Ujj and 
the relative kinetic energy Ty. We calculate the propa- 
gator with the simplified AV8' interaction, the difference 
between the full Hamiltonian H and the simplified H' is 
treated perturbatively. @ 

Because of the large number of spin-isospin amplitudes 
which must be evaluated, the overall speed of the compu- 
tation for larger nuclei is dominated by the calculation of 
the NNN propagator. This term was originally treated 
in a very straightforward manner. Including NNN inter- 
actions, the full propagator G can be calculated as: 



5 



At, 



G^(R,R';At) = exp^r] G (R,R) exp[- ^(^f fc (R) + ^«(R))— ] 



x (a|J 8 (R;^)|7>(7| 



flto,«( r «> r «) 



(14) 



where 



7 3 (R; Ar/2) 



(15) 



Here V^ k is the spin-isospin independent part of the 
NNN interaction, which can be trivially exponentiated. 
The V^ k contains the two-pion-exchange component of 
the NNN interaction, and in the case of the Illinois po- 
tentials, the three-pion-exchange component. 

We have previously employed the fact that the largest 
part of V£ k is a sum of products of spin and isospin 
anticommutators that can be expressed as a sum of terms 



each of which contain only two-nucleon spin and isospin 
operators. This part, which we denote Vijk, is not much 
more difficult to treat than the typical NN potential, i.e., 
while it involves spatial coordinates of three particles, 
the spin-isospin algebra is equivalent to that of an NN 
interaction. 

The remaining terms depend upon the spins and 
isospins of three nucleons and are comparatively weak. 
If we define a propagator G which is obtained from the 
NN potential tty and the Vijk through Eq. (|14|l , the com- 
plete propagator including the full NNN interaction can 
be written as: 



G ai /s(Rw,Ri;MAr) = [1-(%(Rm)- K ife (R M ))M^] 



[qgcr 



At, 



m,Rm-i;At)...G(R 2 ,R i; At) [l-(Vkk(Ri)-V yJk (Ri))M— ], (16) 



where the number M of intermediate steps with the sim- 
plified propagator G is typically four to five. This treat- 
ment is equivalent to evaluating the propagation with 
a simplified Hamiltonian containing only two-body spin 
and isospin operators, and then correcting the end points 
for the difference between the full NNN interaction and 
the simplified one. Clearly this method has the same 
short-time limit as the original implementation, but the 
computational time is significantly decreased because the 
number of full three-nucleon spin-isospin operations is 
much reduced. 

Further generalizing this method, it is possible to in- 
troduce artificial fluctuations into the propagator which 
average to the correct propagator in the short-time limit, 
but are much more efficient computationally. We can re- 
place the factor h{R) above (Eq. IT5|l by 




where 6(x < 0) = and 6{x > 0) = 1. The function 
Pijk is an arbitrary function of the particle coordinates 
(rj,rj,rfe) subject to the condition < Pijk < 1. Inte- 



grating over the "auxiliary fields" trivially recovers 
the original term in the propagator, as with probability 
Pijk the triplet term is Vijk/Pijk and with probability 
1 — Pijk the triplet contribution is zero. The original for- 
mulation for the propagator, Eq. I|15f) . is equivalent to 
choosing P^k = 1. 

Choosing a Pijk which decreases as the three-nucleon 
separation increases is computationally very advanta- 
geous, however. Based upon various trials, we have used 

{1, ^max < Nijk, 

Xijk/Xm&x, -^inin — -^ijk < -^max,(18) 
^min> X-ijk < ^ m j n , 

where 

Xijk=f{r i j)f(ri k )+f{r i j)f(rjk)+f{r ik )f{rjk), (19) 

T(r) = { J M \ r > r ^ (20) 
[_ r(m w r ma x), r < r max . 

Here T{x) is the radial tensor function used in Vijk (sec 
Eqs. (3.5-8) of Ref. 3) and r max is the location of the 
maximum of T{m 1T r). For small triplet separations, as 
effectively specified by X m ax, the full contribution is al- 
ways calculated. For larger separations, the contribution 
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is calculated less often. Because of the l/P^fc in Eq. (|17fl . matrix elements can be computed as: 
one does not want P^ to become too small; X m \ n is 



ijk 

used to control this. Typical values are X max = 5.0 and 

x min = o.oi. 

This technique essentially invokes a larger time step 
for well-separated triplets. Tests are of course required 
to make sure that the overall time step is small enough 
to avoid physically significant errors in the expectation 
values of interest. This method is particularly important 
in nuclear physics applications, where the triplet compu- 
tations are so expensive, but could also be employed in 
other systems where quantum Monte Carlo methods are 
used. 

We have also adopted a similar scheme in computing 
the quadratic spin-orbit (L ■ S) 2 and angular momen- 
tum (L 2 ) contributions to the energy in both VMC and 
GFMC calculations. These terms in the interaction are 
typically quite small at large pair separations, and so we 
can evaluate them with a probability and then mul- 
tiply the calculated contribution with the inverse of Pij . 
Using this scheme enables one to calculate the energies 
and expectation values in the AV18 interaction at a com- 
putational cost much closer to evaluating the simplified 
AV8' interaction. 



Nh{t) 



[HijHji] 2 , 



(24) 
(25) 



where we use expectation values computed from sepa- 
rately propagated *j(r) and ^(r). For i = j these 
equations reduce to the standard GFMC calculation de- 
scribed in Sec. III. 

For larger nuclei, A > 8, the fermion sign problem, i.e., 
the problem of symmetric noise, becomes a sufficiently 
large computational burden^hat a "constrained path" 
algorithm must be utilized In this case the propa- 
gator for different states is not identical, as it involves 
constraints based upon the different trial wave functions, 
and the above equations are only approximate. As is de- 
scribed in 3] for calculations of the lowest state of given 
(J^, T), any errors in the energy that are introduced by 
the constrained path can be removed by doing 10 to 20 
unconstrained steps before evaluating the energy. In the 
present work we follow this same proceedure and make 
10 to 20 unconstrained steps before evaluating Hij and 

Ni4. 



IV. GFMC EVALUATION OF EXCITED 
STATES 



It is possible to treat at least a few excited states 
with the same quantum numbers using VMC and GFMC 
methods. The VMC calculations have been described 
above, and essentially involve solving a generalized eigen- 
value problem, Eqs. © and J7J). The same basic method 
can be applied in GFMC calculations, though the im- 
plementation is slightly more involved. In this section, 
^T.i represents the trial wave function for the i th state 
of specified (J n ,T) and ^(r) is the GFMC wave function 
propagated from it. By construction (^taI^tj) = for 
i 7^ j. We would like to calculate the Hamiltonian and 
normalization matrix elements as a function of r: 



Hijir) 



(*Mt/2)| H |^-(r/2)) 
1*^/2)11*^/2)1 : 
(^(7-/2)1^(7-/2)) 
|* < (t/2)||* 3 .(t/2)|' 



(21) 
(22) 



where |*j 



|('I'i|*i)| 1 / 2 . Solving the generalized eigen- 
value problem with these Hamiltonian and normalization 
matrix elements would yield improved upper bounds for 
the ground and low-lying excited states of the system. In 
the limit r — > oo the solutions would be exact. 

In GFMC we can compute mixed expectation values 
such as 



(* t (r)| Q \*rj) 
(*i(r)|*T,i) : 



(23) 



where the denominator involves just state i. Since the 
propagator commutes with the Hamiltonian, the desired 



-15, 
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7 Li(5/2") 
AV18 + Ifflnois-2 
GFMC Propagation 



^ 4*5/2" 

♦ 3 rd 5/2- 
■ 2 nd 5/2" 

• I st 5/2" 




0.4 0.6 

T (MeV" 1 ) 

FIG. 1: (Color online) GFMC energies of four § states in 
7 Li versus imaginary time, r. The solid symbols show the 
computed energies at each r; open symbols show the results 
of the rediagonalization discussed in the text. 

As an example, Figs. ^ and [21 show the computation 
of the energies of four | states in 7 Li; the *T,i used to 
start the propagations are defined in Table Hvl The di- 
agonal Ha (t) are shown as solid symbols in Fig. ^ The 
lowest state has mainly [43] symmetry and can easily de- 
cay to the a+t channel; it has a large experimental width 
(918 keV) and its computed energy is slowly decreasing 
to the energy of the separated clusters. The remaining 
states are mainly of [421] symmetry and so are princi- 
pally connected to the 6 Li+n channel. The second | 
state is experimentally just above the 6 Li+n threshold 
and has a small width (80 keV); its computed energy 
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0.2 0.4 0.6 0.8 



X (MeV- 1 ) 

FIG. 2: (Color online) Overlaps of the GFMC wave functions 
for the first | state in 7 Li with the other three | states 
versus imaginary time. 

becomes constant with increasing r. The last two § 
states are not experimentally known, but the very slow 
decrease with r of the energy of the third state suggests 
that this state may also be narrow. Figure [21 shows the 
off-diagonal overlaps Nh(t); they are small and do not 
show signs of steadily increasing with increasing r. The 
Ni-tj (t = 0) are not identically zero because the diagonal- 
ization that determined the Plsw was made in a differ- 
ent Monte Carlo walk from the ones generating the Nij . 
These results show that the (constrained) GFMC prop- 
agation largely retains the orthogonality of the starting 
^T,i- Contrary to what might have been expected, the 
propagation of the higher states does not quickly collapse 
to the lowest state. 

The open symbols in Fig. ^ show the results of gener- 
alized eigenvalue calculations made for each r using the 
Hij(r) and ATy(r). As is expected from the small over- 
laps, there is no statistically significant change of the re- 
diagonalized energies from the directly computed Hu{t). 
Based on this result and the generally constant Hu(t) 
that we obtain for other cases, we use the Hu(r) directly 
without making a final rediagonalization. 



TABLE VIII: Energies and excitation energies (MeV) of se- 
lected nuclear states calculated for the AV18/UIX Hamilto- 
nian by VMC and GFMC. Monte Carlo statistical errors in 
the last digits are shown in parentheses. See the text for 
italicized GFMC values. 



Energy Excitation Energy 





VMC 


GFMC 


VMC 


GFMC 


*He(0 + ) 


-27.39(3) 


-28.31(2) 






6 He(0+) 
6 He(2+) 


-24.79(9) 


-28.02(9) 






-23.07(9) 


-26.08(9) 


1.7(1) 


1.9(1) 


6 He(l+) 


-20.72(9) 


-23.9(2) 


4.1(1) 


4.2(2) 


6 Li(l+) 


-27.96(8) 


-31.25(8) 






6 Li(3+) 


-25.05(9) 


-28.45(8) 


2.9(1) 


2.8(1) 


6 Li(2+) 


-23.8(1) 


-27.25(8) 


4.2(1) 


4.0(1) 


6 Li(l + , 2 nd ) 


-22.69(8) 


-26.2(1) 


5.3(1) 


5.1(1) 


7 He(|-) 


-20.93(8) 


-26.3(1) 






7 He(i") 


-18.8(1) 


-25.2(1) 


2.1(1) 


1.1(1) 


7 He(D 


-18.33(9) 


-23.9(1) 


2.6(1) 


2.5(2) 


7 Li(§~) 


-33.0(3) 


-37.5(1) 








-32.9(1) 


-37.6(1) 


.1(3) 


-.1(2) 


7 Li(|-) 


-27.2(1) 


-32.2(1) 


5.8(3) 


5.3(1) 


7Li (n 


-26.61(9) 


-31.1(1) 


6.4(3) 


6.5(2) 




—23.7(1) 


—29.7(2) 


9.3(3) 


7 8(2) 




— 99 Sfl 1 


—29.1(2) 


1 n 9(X\ 

l-\J.z.yo J 


o.oyoj 


Mo" ' 2 ) 


01 ( 1 ^ 

— Zl.8(l ) 


—ZQ.L{Z ) 


LL.Z(o) 


y.o{z ) 


7 TU 5 ~ 1 rd \ 
Ll (2 ' 6 ) 


on 


97 f\(0\ 

— z i .y}\Z) 




in Kfo\ 


7 Tif 5_ d th \ 
^Ka ' 4 1 




9/ Id} 




1 . 1 [U J 


8 He(0+) 


-20.75(7) 


-27.7(1) 






8 He(2+) 


-18.2(1) 


-25.0(1) 


2.5(1) 


2.8(2) 


8 He(l+) 


-17.56(8) 


-23.3(3) 


3.2(1) 


4-4(3) 


8 Li(2+) 


-30.7(1) 


-38.8(1) 






8 Li(l+) 


-29.6(1) 


-37.8(2) 


1.1(1) 


1.0(2) 


8 Li(0+) 


-28.6(1) 


-36.9(1) 


2.1(1) 


2.0(2) 


8 Li(3+) 


-27.24(9) 


-35.4(1) 


3.5(1) 


3.4(2) 


8 Li(l+, 2 nd ) 


-27.7(1) 


-35.2(1) 


3.0(2) 


3.6(2) 


8 Li(4+) 


-24.42(9) 


-32.3(2) 


6.3(1) 


6.5(2) 


8 Be(0+) 


-48.6(1) 


-55.2(1) 






8 Be(2+) 


-45.71(9) 


-52.1(2) 


2.9(1) 


3.1(3) 


8 Be(4+) 


-37.59(9) 


-412(5) 


11.0(1) 


11.0(5) 


8 Be(l+) 


-28.07(8) 


-37.0(2) 


20.6(1) 


18.2(2) 


8 Be(2+, 2 nd ) 


-28.8(1) 


-36.4(6) 


19.8(1) 


18.9(6) 


8 Be(3+) 


-26.69(7) 


-35.2(2) 


21.9(1) 


20.0(2) 



V. NUMERICAL RESULTS FOR A = 6 TO 8 
NUCLEI 

Figure^shows that for some unbound cases the GFMC 
energy never stabilizes, but rather slowly decays with in- 
creasing r. In these cases a linear fit with non-zero slope 
provides a much lower \ 2 fit to the E(r) for large r than 
does a constant fit. We use this linear fit to extrapolate 
back toT« 0.08 because it appears that the GFMC cal- 
culations for bound systems have all become stable by 
t w 0.08. Such extrapolated values are printed in ital- 
ics in the tables; they must be considered as less reliable 
than the non-italicized values. 

Selected VMC and GFMC energies computed for the 



AV18/ULX Hamiltonian are shown in Table IVnTl The 
VMC energies for the ^>t of Eq. (JTJ are significantly im- 
proved compared to those reported in Refs. The 
4 He energy has been lowered by 0.4 MeV, the A=Q en- 
ergies by 1 MeV, the A=7 energies by 2 MeV, and the 
A = 8 energies by 2.75 MeV. Compared to the more 
sophisticated (and more expensive) 'Fy of Refs. 0, El, 
which includes spin-orbit and additional three-body cor- 
relations, the energies from the present 'Ft are as good 
for A— 6 and even better for A— 7,8. The improvements 
are due primarily to 1) improved shapes for the two-body 
tensor correlation functions ut(r), u tT (r) contained in f/y 
[Eq. J3J], including allowing them to be different for s- 
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TABLE IX: A = 4, 6 and 7 GFMC energies for the AV8', AV18, and AV18/IL2 Hamiltonians compared with experimental val- 
ues [Til ]. Except as noted, states have the same isospin as the corresponding ground state. Square brackets enclose experimental 
energies which have uncertain J 71 assignments. See the text for italicized GFMC values. All values are in MeV. 



AV8' AV18 AV18/IL2 Experiment 

Energy Width 



*He(0 + ) 


-25.14(2) 


-24.07(4) 


-28.37(3) 


-28.30(0) 




6 He(0 + ) 


-25.11(3) 


-23.8(1) 


-29.28(2) 


-29.27 




6 He(2+) 


-23.21(4) 


-21.9(1) 


-27.3(2) 


-27.47(3) 


0.113 


6 He(2+, 2 nd ) 


-81.4(1) 


-20.4(1) 


-24.6(1) 






6 He(l + ) 


-21.0(1) 


-19.6(2) 


-24.5(1) 






6 He(0 + , 2 nd ) 


-20.15(6) 


-19.0(1) 


-23.3(1) 






6 Li(l+) 


-28.15(3) 


-26.9(1) 


-32.0(1) 


-31.99 




6 Li(3+) 


-25.33(3) 


-24.1(1) 


-29.8(1) 


-29.80 


0.024 




0/1 Q1 (1\ 

— z4.oi^oJ 


— Zo.\J[i ) 




OQ A Q 

— zo.4o 




6 Li(2+) 


-24.18(5) 


-22.8(1) 


-27.8(1) 


-27.68(2) 


1.3 


6 Li(2+, r=i) 


-22.48(4) 


-21.1(1) 


-26.5(1) 


-26.62 


0.54 


6 Li(l+, 2 nd ) 


-23.19(6) 


-22.0(1) 


-26.4(2) 


-26.34(5) 


1.5 


6 Li(l+, 3 rd ) 


-19.7(2) 


-18.8(1) 


-23.4(1) 






7 He(I~) 


-23.39(7) 


-21.9(1) 


-28.8(2) 


-28.82(3) 


0.16 


7 He(±~) 


-22.05(9) 


-20.7(1) 


-25.9(2) 






7 He(|~) 


-20.8(1) 


-19 1(1 ) 


-25.5(1) 


[-25.90(10)] 


1.99 


7 He(^~ 2 nd ) 






—25.0(1) 






7 Hef£~ 3 rd ) 






-21.7(3) 






Ll l2 > 


oo oat a \ 
— oo. 00^4 J 


qo C\( 1 ^ 
— OZ.UI^ 1 ) 


— oo.y^i ) 


QO> O/l 

— oy.z4 




7 Li(I-) 


-33.95(4) 


-32.2(1) 


-38.7(1) 


-38.76 




7 Li(l-) 


-28.6(1) 


-26.8(2) 


-34.0(1) 


-34.59(1) 


0.069 


7 Li(D 


-28.11(6) 


-26.4(1) 


-32.3(1) 


-32.64(5) 


0.92 


7 Li(f ~, 2 nd ) 


-26.42(7) 


-24.5(1) 


-31.7(1) 


-31.79 


0.08 


7 Li(|", 2 nd ) 


-25.4(2) 


-24.4(1) 


-29.7(1) 


-30.49 


4.7 


7 Li(±~, 2 nd ) 


-25.37(7) 


-23.8(1) 


-29.1(2) 


-30.15 


2.8 


7 Li(|", 2 nd ) 


-24.59(7) 


-23.0(1) 


-29.2(2) 


-29.67(10) 


0.44 


7 Li(§~, 3 rd ) 






-28.6(1) 






7 Li(| _ , 4 th ) 






-25.8(1) 







shell and p-shell nuclei, 2) letting the f S p{r) vary with 
the LSfri] wave function component, and 3) using more 
extended f S p(j") and f PP {r) which allow the overall wave 
function to be more diffuse. In particular, the lack of 
experimental charge radii for A=8 nuclei makes it diffi- 
cult to fix the optimal size of these wave functions; in 
the present case the matter radii for the ground states 
have been adjusted to match the earlier GFMC results 
of Ref. Q. 

Many of the GFMC results in Tablc lVlTll havc changed 
by 1 to 3 standard deviations from the results reported 
in Ref. Q|. However, the quoted errors are statistical 
only; we have previously estimated 0, 01 the systematic 
errors to be of the order of 1-2%, and the changes re- 
ported here are within that range. The changes are due 
to the improved treatment of Vijk discussed above, prop- 
agation to larger r, and in some cases improved Plsw 
choices. These changes are significantly less than the cor- 
responding improvements in the \I/t values, indicating 
that (aside from modifications of the Plsm) high-energy 
contaminations have been removed from the \&t; these 
were previously easily corrected by the GFMC. 



For the p-shell nuclei, the improved energies are 
still 10 to 20% higher than the GFMC values. How- 
ever the excitation energies computed with the are 
generally quite accurate, providing the states are domi- 
nated by the same spatial symmetry as the ground state. 
The second | and higher states of 7 Li and the 1 + and 
3 + states of 8 Be have lower spatial symmetry than the 
corresponding ground states and their VMC excitation 
energies are significantly too high; we also see this pat- 
tern in the VMC energies of other higher states and with 
other Hamiltonians. 

Tables |IX] and H give GFMC energies for the AV8', 
AV18, and AV18/IL2 Hamiltonians and also show the 
experimental energies and widths of the states that 
have been observed (widths less than 0.01 MeV are not 
shown). Figure |3 shows the corresponding excitation en- 
ergies. The values for second or higher states of a given 
J n are presented here for the first time, however the other 
values have all been recomputed with, in many cases, 
increased maximum t, improved determinations of the 
(3ls[ti] in the ^t, and improved treatment of the Vijk 
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TABLE X: A = 8 GFMC energies for the AV8', AV18, and AV18/IL2 Hamiltonians compared with experimental values |Tl|. 
Conventions are the same as in Table HX| The 2+ states near —39.7 MeV in 8 Be are strongly isospin mixed; their experimental 
energies, marked with an asterisk, have not been corrected for this effect. 





AV8' 


AV18 


AV18/IL2 


Experiment 
Energy Width 


s He(0 + ) 


-24.3(1) 


-23.0(1) 


-31.72(4) 


-31.41(1) 




8 He(2+) 


-21.93(6) 


-20.4 1) 


-27.0(2) 


-28.31(50) 


0.6 


8 He(l+) 


-20.3(2) 


-18.8(3) 


-25.9(2) 






8 He(0 + , 2 nd ) 






-24.8(2) 






8 He(2+, 2 nd ) 






-23.7(2) 






Ll(2 ) 


OA *7A((\\ 

— o4. 1 4(0 ) 


00 n( 1 \ 
— 6Z. / ^ ± J 


/II Q/0\ 


-41.28 






—34.34(7) 


32.1(1) 


—40.5(2) 


-40.30 




8 Li(3+) 


-32.16(7) 


-30.1(2) 


-39.4(2) 


-39.02 


0.033 


8 Li(0+) 


-33.51(6) 


-31.5(1) 


-38.3(2) 






8 Li(l+, 2 nd ) 


-32.76(7) 


-31.0(1) 


-37.2(2) 


-38.07 


1. 


8 Li(2+, 2 nd ) 


-31.2(2) 


-29.7(1) 


-36.6(4) 






8 Li(2+, 3 rd ) 


-30.49(7) 


-28.7(2) 


-36.9(2) 






8 Li(l + 3 rd ) 


-31.14(7) 


-29.1(2) 


-35.9(2) 


[-35.88] 


0.65 


8 Li(3+, 2 nd ) 


-28.47(7) 


-26.3(1) 


-33.8(2) 


[-35.18(10)] 


1. 


8 Li(4+) 


-29.1(1) 


-27.1(1) 


-34.7(2) 


-34.75(2) 


0.035 




AO n CT / T\ 

—48.95(7) 


—46.3(2) 


—56.3(1) 


—56.50 




Be(2 T ) 


—45.6(3) 


—43.7(2) 


— 52.2(2) 


CO A A 1 0\ 

— 00. 44(3) 


1.4 


Be(4 J 


— of. 1(0) 


or if'o'* 
— oO.Z{Z ) 


AZ. A { 0\ 

— 40.4( ) 


— 40.10(10 ) 


. 


8 "p / Y) + rp i\ 

r>e^z , 1 — 1 J 


QQ Af\( a\ 


Q1 0/1 

— ol.z(l J 


/in a (o\ 

— 4U.4(Z ) 


QQ 07* 

— oy.o I 


U.lUo 


raeyZ , Z ) 


QQ A 1 ((\\ 

— oo.41(y ) 


— ol.U(z J 


— oy.o(o ) 


oy.oo 


U.U I 4: 


8D./1+ rp i \ 

r>e^i , i — i j 


QQ f\of 7^ 


on 7( -\\ 


— Oif.UyZ ) 


-38.86 


0.011 


8 Be(l+) 


— ^ 9^(q\ 


30.8 (2) 


—39.4(3) 


-38.35 


0.138 


8 Be(3+, T=l) 


-30.87(7) 


-28.8(2) 


-38.0(2) 


-37.43 


0.27 


8 Be(l + , 2 nd ) 






-37.3(3) 






8 Be(3+) 


-31.68(9) 


-29.7(2) 


-37.2(1) 


-37.26(1) 


0.23 


8 Be(4+, 2 nd ) 






-37.6(3) 


-36.64(5) 


0.7 


8 Be(2+, 3 rd ) 


-31.49(9) 


-29.1(2) 


-37.0(3) 


-36.40 


0.88 


8 Be(0 + , 2 nd ) 






-36.0(3) 


-36.30 


0.72 


8 Be(3+, 2 nd ) 






-34.0(3) 


-35.00 


1. 


8 Be(2+, 4 th ) 






-34.0(3) 


-34.30 


0.8 



in the propagation. For these reasons, some of the val- 
ues are significantly different from previously published 
values. 

The AV18/IL2 Hamiltonian was originally determined 
by making a three-parameter fit to 17 states in the A = 3 
to 8 region [4(. There are 36 states with 6 < A < 8 that 
have experimental energies in Tables llXl andlXl The RMS 
error in predicting the energies of these states is only 0.60 
MeV; the RMS error for just the 17 states with widths 
less than 0.2 MeV is 0.38 MeV and the RMS error for the 
7 ground-state energies is 0.31 MeV. The RMS errors 
in excitation energies are 0.76 MeV for all 29 excited 
states and 0.54 MeV for the 10 narrow states. Analog 
states have been omitted from all of these averages but 
the states with uncertain experimental assignments 
were included; omitting the latter does not substantially 
change the above numbers. These results are basically 
the same as the RMS errors reported for 26 states in 
Rcf. 0. 

Extensive breakdowns of the total energies were pre- 
sented in Ref . . The improved treatment of V«fc in the 
GFMC propagation has resulted in significant changes 



TABLE XI: 


GFMC expectation values 


of Vij 


AV18/IL2 Hamiltonian 










y37r yR 


Vijk 


4 He 


-16.3(1) 


0.63(1) 7.3(1) 


-8.4(1) 


6 He 


-18.9(1) 


-0.58(1) 8.9(1) - 


-10.6(1) 


6 Li 


-20.1(2) 


-0.09(3) 9.4(2) - 


-10.8(2) 


7 He 


-23.5(3) 


-2.22(6) 11.6(2) - 


-14.2(2) 


7 Li 


-24.9(2) 


-0.50(3) 11.7(1) - 


-13.7(2) 


8 He 


-27.0(1) 


-4.36(2) 13.8(1) - 


-17.5(1) 


8 Li 


-31.4(3) 


-2.6(1) 15.8(2) - 


-18.2(3) 


8 Be 


-35.9(2) 


0.32(3) 16.4(2) - 


-19.2(2) 



to the expectation values of Vijk and its components for 
the AV18/IL2 Hamiltonian; revised values are shown in 
Table IXII The other contributions to the total ener- 
gies shown in Ref. £j are not significantly changed. The 
isovector and isotensor energies presented there have also 
not significantly changed in the current calculations. Ta- 
bles [DC] an d K] contain energies of a few T=l states for 
6 Li and 8 Be; these have been computed perturbatively 



10 



25 



20 



15 



> 



P9 10 



CD 

d 







Argonne v 18 
With Illinois-2 




5/2" 






+ - 



xs 

re 



Op' 



FIG. 3: (Color online) GFMC excitation energies computed for the AV18, AV18/UIX, and AV18/IL2 Hamiltonians compared 
with experimental values. The shaded bands show the Monte Carlo statistical (or experimental) errors. The narrow shaded 
bars on the experimental values show the experimental widths. 



using the isovector and isotensor energies computed for 
6 He and 8 Li, respectively. Recently we have realized 
that there is a significant difficulty in extracting pre- 
cise RMS radii from GFMC calculations, especially for 
weakly bound systems. For this reason we are deferring 
presenting updated values for the RMS radii. 



nucleon force in the GFMC calculations. This has altered 
the detailed breakdown of NNN contributions as shown 
in Tablc lXll We expect a similar alteration of these terms 
in the A=9,10 nuclei, but the total energies should stay 
within 1-2% of our previously reported values |5(. 



VI. CONCLUSIONS 
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to compute the energies of multiple nuclear states with 
the same quantum numbers. This substantially increases 
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